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1.  Introduction 


Future  military  analog  fiber-optic  links  require  high-bandwidth  (20-30  GHz)  and  large 
dynamic  range  (>120  dB)  that  are,  furthermore,  compact,  light  weight,  and  with  low 
power  consumption.  Technologies  that  have  been  investigated  so  far  can  not  meet  the 
above-mentioned  demands.  These  technologies  include,  cw  laser  diode  direct 
modulation,  mode-locked  laser  diode  external  modulation,  and  external  modulation. 

Optical  coupling  of  a  Laser  diode  (LD)  with  an  external  cavity  (EC),  such  as 
shown  schematically  in  Fig.  1,  has  been  used  to  improve  laser  frequency  stability  and  to 
narrow  laser  linewidth.  Such  coupling  can  also  significantly  change  LD  dynamics.  For 
long  ( >  1  cm)  external  cavity  LDs  these  changes  lead  to  chaotic  self-modulation  of  laser 
light  with  dramatic  broadening  of  laser  spectrum  ("coherence-collapse"  effect).  For  a 
short  EC  LD  there  exists  another  type  of  instability  associated  with  the  EC  mode 
interactions  alone  [1].  The  compound-cavity  mode  interaction  inside  the  regime  of 
stable  CW  operation  lead  to  the  generation  of  high  frequency  (>20  GHz)  microwave 
oscillations.  Such  high-frequency  oscillations  arise,  as  has  been  conjectured  [1],  if  the 
cavity  losses  for  the  main  and  side  compound-cavity  modes  of  the  short  EC  LD  become 
equal  with  increasing  feedback.  This  equalization  does  not  lead,  however,  to  a  switching 
of  lasing  to  this  sidemode,  but  instead  induces  a  self-mode-locking  regime  with  deep 
oscillations  of  the  laser  output. 
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We  propose  to  use  short  EC  LDs  as  tunable  microwave  oscillators  with  frequency 
modulation  capabilities.  There  are  two  way  the  frequency  can  be  modulated  [2] . 

•  It  has  been  shown  that  the  oscillation  frequency  increases  with  increasing 
feedback.  Therefore,  the  frequency  of  the  microwave  oscillation  can  be 
modulated  by  the  feedback 

•  It  should  be  noted  that  the  oscillation  is  a  result  of  the  interaction  of  compound- 
cavity  mode  interaction.  Therefore,  once  the  optical  length  of  the  individual 
cavity  is  changed,  the  oscillation  is  tuned. 

We  have  previously  computed  the  dynamics  of  short  EC  LDs,  mapped  out  the  system 
parameter  space  for  microwave  oscillations,  and  determined  the  frequency  range  of  the 
microwave  oscillations. 

The  ultimate  research  objective  is  to  experimentally  demonstrate  the  generation  of 
tunable  microwave  intensity  oscillations  in  laser  diodes  with  short  external  cavity  for  use 
in  analog  fiber-optic  links.  A  better  understanding  of  the  fundamental  mechanism  is 
therefore  essential  for  progress  towards  employing  these  devices  in  real  applications.  In 
this  report,  we  document  enhanced  understanding  of  the  underlying  physical  principle 
for  the  microwave-frequency  intensity  modulation  from  a  nonlinear  dynamics  point  of 
view,  and  ways  to  experimentally  measure  the  noise  characteristics  of  the  oscillations. 
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2.  Formulation  of  the  Problem 


2.1  Rate  Equation 

Laser  Diode  cavity  External  Cavity 

Rl  R2  R3 


p,  _  ■ 

J 

| Kliiiii  S  111 

fextP  1 

L _ _ L _ _ J 

r  i 

L 

l 

Lext 

Fig.  1  Schematic  of  a  laser  diode  with  an  external  cavity. 

The  schematic  of  a  semiconductor  laser  diode  with  an  external  cavity  is  shown  in 
Fig.  1.  Rl  and  R 2  are  the  reflectivity  of  the  laser  facets  and  R3  is  the  reflectivity  of  the 
external  mirror.  L  and  Lext  are  the  lengths  of  the  laser  diode  cavity  and  the  EC;  fext  is 
the  feedback.  We  will  consider  tcd*  <<1  ,  where  r=  2  Le„  /  c  is  the  photon  round-trip 
time  in  EC  and  coR  is  the  angular  frequency  of  the  relaxation  oscillation. 

We  adopt  the  Lang  and  Kobayashi  rate  equations  for  the  SEC  LD.  There  are 
several  assumptions  for  these  rate  equations.  First,  the  'solitary  laser  is  assumed  to  be 
single  mode,  which  implies  r»  r0  (r0  =  2Lng  I  c  ,  is  the  photon  round-trip  time  in  the 

LD-cavity,  ng  is  the  group  refraction  index  of  the  LD.).  Second,  fexl  «  1  which  means 
only  one  round-trip  of  the  emitted  photons  is  significant  inside  the  EC.  Then  we  arrive  at 
the  following  two  rate  equations  : 
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m) 

dt 


=  J 


ns-Q+ZLgN(N-Ns) 


E(t)  +  Ug~—)E(t )  +  KE(t  -  r)exp(-zQr) 

1  *D  U) 


dN  N 

—  =  J - gS 

dt  t 


(2) 


where  g  =  g(N)(\ -  KpP)  is  the  nonlinear  gain  where  the  differential  gain  is  given  by 
g(N)  =  gN  (N  -  Ns),  P  is  the  output  power,  Qf  is  the  resonance  optical  frequency  of 
the  LD  without  EC,  a  is  the  linewidth  enhancement  factor,  gN  is  the  differential  gain, 
tp  is  the  photon  lifetime  in  the  LD  cavity,  r  is  the  carrier  lifetime  due  to  spontaneous 

and  nonradiative  recombination,  5'  =  |£(0|2  is  the  photon  number  in  the  LD  cavity, 
(output  power  P~S),  J  is  the  carrier  injection  rate.  The  feedback  parameter  is  given  by 

C, 


K  = 


where  parameter  C,  which  is  a  function  of  the  facet  and  external  mirror  reflectivity  R2 
and  R2 

C,=(l-*2) 

is  a  measure  of  the  coupling  strength  between  the  two  cavities.  The  rate  equations  for 
the  EC  LD  are  given  in  equations  (1)  and  (2).  The  absolute  squared  value  of  the  electric 
field  amplitude  corresponds  to  the  number  of  photons  in  the  laser  cavity.  Equations  (1) 
and  (2)  can  be  cast  into  rate  equations  for  the  photon  number  S,  the  phase  (f>  ,  and  the 
carrier  number  N: 
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2.2  Numerical  Simulations 

In  numerical  simulations,  the  system  of  rate  equations  system  (3)-(5)  is  solved  by 
using  a  fifth-order  Runge-Kutta  method  [3], [4].  We  specify  the  initial  conditions 
S(0)  =  10,  0(O)=O  and  N(0)=0.  In  the  time  interval  0<t<  r,  the  data  of  photon 
number  S  and  phase  <f>  are  stored  in  memory  and  they  are  reused  after  the  round-trip 
time  t  of  the  EC  in  order  to  account  for  the  optical  feedback. 


0.00  500.00  1000.00  1500.00 


T  im  e  (p  s  ) 


Fig.  2  The  strongest  oscillation  for  the  Lext=5.25mm  and  the  fexl  group.  The 
feedback  parameter  is  0.206.  The  oscillation  frequency  is  equal  to  21.3  GHz. 
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2.3  An  Example  of  High-Frequency  Microwave  Oscillation 


For  example,  we  set  the  output  power  P =5  mW/per  facet.  The  calculations  used  the 
laser  parameters  given  in  Table  I.  Here,  we  also  choose  the  EC  length  Lexl  = 5.25  mm. 
We  found  that  the  feedback  parameters  for  the  strongest  modulation  index  occurs  when 
fext  =  0.206.  Further  increasing  fex,  destroys  the  twin  mode  locking.  It  leads  the  laser 
operation  on  a  new  single-mode  state  with  higher  output  power  (7.06  mW).  The 
strongest  oscillation  frequency  we  get  in  the  group  fex,  is  21.3  GHz  as  shown  in  Fig.  2. 


TABLE  I 

List  of  Parameters  for  Fabry-Perot  laser 


Parameter 

Value 

Units 

scattering  loss  in  the  active  region 

3.6  xlO3 

m~l 

carrier  number  at  transparency 

2.0  x 108 

output  power 

P" 

2.906  x  10~8  x  S 

mW  /  facet 

index  of  active  region 

ns 

4.5 

differential  gain 

Sn 

2.2  xlO3 

s -> 

linewidth  enhancement  factor 

a 

4 

laser  wavelength 

X 

1.3 

um 

spontaneous  emission  lifetime 

T 

2.0  xlO'9 

s 

photon  lifetime 

mm 

1.29072  x  10"12 

s 

gain  saturation  coefficient 

K 

4.8 

w-' 

3.  Bifurcation  Analysis  of  the  ECLD 


Since  the  problem  is  nonlinear  in  nature,  there  does  not  exist  a  general  analytical 
methodology  that  can  lead  to  a  global  conclusion  of  the  operation  characteristics  of  the 
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ECLD.  Our  approach  is  to  establish  a  numerical  technique  that  is  capable  of  predicting 
the  nature  of  the  ECLD  output  given  a  complete  set  of  parameters  for  the  ECLD 
without  numerical  integration  described  in  Section  2.  We  used  a  public-domain  software 
called  AUTO  to  investigate  the  dynamics  of  laser  diodes  with  short  external  cavity. 
AUTO  was  developed  in  the  Applied  Mathematics  Department  at  the  California 
Institute  of  Technology  as  one  of  the  projects  of  the  Differential  Equations  Group  in  the 
Center  for  Research  on  Parallel  Computation.  In  what  follows,  we  describe  the  interface 
FORTRAN  code  for  AUTO  to  simulate  the  laser  diodes  with  external  cavity  and  some 
representative  results.  We  stress  that  the  results  shown  here  only  apply  to  the  particular 
parameters  chosen  as  indicated.  To  find  the  dynamic  behavior  of  ECLDs  with  different 
parameters,  one  needs  to  run  the  AUTO  program  with  the  interface  FORTRAN  code 
with  the  modified  ECLD  parameters. 

3.1  Bifurcation  Diagrams  of  ECLD  Steady  States 

To  construct  the  bifurcation  diagrams  of  ECLD  steady  states,  the  first  step  is  to  find  all 
ECLD  steady  states  for  a  given  set  of  the  ECLD  parameters  and  then  using  AUTO  to 
finding  the  continuations  of  the  steady  states. 

To  find  all  ECLD  steady  states  for  a  given  set  of  the  ECLD  parameters,  one 
needs  to  find  all  the  solutions  to  the  nonlinear  algebraic  equations: 
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(6) 


(g  -  —)S(  t)  +  2^V^(0V^-7)  C0S(Q, T  +  0(0  -  0('  -  0)  =  o 

rP 

J-  — -g5(0  =  0  (7) 

Tsp 

arrived  by  setting  the  right-hand  side  of  Eqs.  (3)  and  (5)  to  be  equal  to  zero. 

Once  all  the  steady  states  of  the  ECLD  for  the  given  set  of  parameters  are  found, 
the  bifurcation  diagram  can  be  constructed  by  continuing  the  steady  states  as  one  of  the 
parameters  in  the  parameter  set  is  varied  in  both  directions  (increasing  or  decreasing).  It 
is  found  that  Hopf  bifurcation  delineates  stable  and  unstable  steady  states. 

3.2  Stability  Boundaries  of  ECLD  Modes 

To  construct  the  Stability-Boundary  diagram  of  ECLD  steady  states,  the  first  step  is  to 
find  the  Hopf  bifurcation  points  of  all  ECLD  steady  states  for  a  given  set  of  the  ECLD 
parameters  and  then  using  AUTO  to  finding  the  continuations  of  the  Hopf  bifurcation 
points. 


To  find  Hopf  bifurcation  points  of  all  ECLD  steady  states  for  a  given  set  of  the 
ECLD  parameters,  one  needs  to  find  all  the  solutions  to  the  nonlinear  algebraic 
equations  (6),  (7),  and  (8) 

D(s)  =  -s3  +  s2  [r  +  2rce(r)  cos  <f>]  -  s[<y  R2  +  ic2e2  (r)  +  2y  eKe(z)  cos  <z )] 
coR2e(T)( cos <f>  -  a  sin  <j>)  +  y  eKe2  (z)  +  ay  j  pe(z)  sin 


+K 


(8) 
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where  e(r)-l-e,<ur  and  <j)-Q r.  Here  ye  =1/  rs  +  gNS  is  the  carrier  perturbation 
decay  rate,  y  =  ye  +  yp  is  the  decay  rate  of  relaxation  oscillations,  with  y p  =  gKpP 
being  the  contribution  of  the  nonlinear  gain.  The  angular  frequency  of  relaxation 
oscillations  coR  is  given  by  a>R  =  ggNS  +  yey p.  Equation  (8)  is  arrived  by  setting  the 
determinant  of  the  linearized  ECLD  system  to  be  equal  to  zero. 

SQRT  [FEXT] 


-0.015  -0.010  -0.005  0.000  0.005  0.010  0.015 

TD-50  [PS] 

Fig.  3  The  stability  boundaries  of  4  ECLD  modes  in  the  -  rd  plane.  The  injection 
current  is  at  twice  the  threshold.  The  meaning  of  LP  is  explained  subsequently. 

Once  Hopf  bifurcation  points  of  all  the  steady  states  of  the  ECLD  for  the  given 
set  of  parameters  are  found,  the  bifurcation  diagram  can  be  constructed  by  continuing 
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the  steady  states  when  one  of  the  parameters  in  the  parameter  set  is  varied  in  both 
directions  (increasing  or  decreasing).  An  example  is  shown  in  Fig.  3. 


3.3  Regions  of  Existence  of  ECLD  Modes 
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Fig.  4  Limiting  points  (LP)  of  existence  of  modes  superimposed  on  the  stability 
boundaries.  When  LP  of  a  higher-order  mode  is  above  the  stability  boundary,  robust 
high-frequency  oscillations  result.  Parameter  same  as  in  Fig.  3,  except  that  injection  is 

three  times  the  threshold. 


It  is  found  by  numerical  simulation  that  even  in  regions  above  the  stability 
boundary  of  a  particular,  the  high-frequency  oscillations  are  not  always  present.  We 
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believe  that  this  due  to  the  existence  of  other  modes,  which  might  be  stable  in  the 
unstable  region  of  the  mode  that  we  are  interested  in.  To  do  so,  we  calculated  the 
regions  of  existence  of  all  modes.  High-frequency  oscillations  are  produced  only  if  the 
boundary  of  existence  [limiting  point  (LP)  of  existence]  of  higher-order  modes  is  above 
the  stability  boundary.  An  example  where  robust  high-frequency  oscillations  are  possible 
is  shown  in  Fig.  4.  A  careful  examination  of  Fig.  3  reveals  that  robust  high-frequency 
oscillations  are  not  possible,  even  though  the  bifurcation  analysis  alone  predicts  the 
existence  of  the  high-frequency  oscillations. 

3.4  Experimental  Generation  of  High-Frequency  Oscillation 

Using  a  Fabry-Perot  laser  diode  and  an  optical  flat  as  the  external  mirror,  The  high- 
frequency  intensity  oscillations  have  been  generated  in  Rome  Lab  by  Steve  Johns.  The 
highest  frequency  achieved  is  9.5  GHz.  Modulation  of  these  oscillations  has  also  been 
demonstrated.  Up  to  750  MHz  modulation  has  shown  excellent  frequency  modulation 
characteristics  with  minimal  residual  amplitude  modulation. 
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4.  Investigation  of  the  Noise  Characteristics 

4.1  Preliminary  Comment 

Periodic  orbits  in  nonlinear  optical  systems,  such  as  the  high-frequency  oscillations 
described  in  Sections  2  and  3  can  serve  as  subcarriers  in  optical  communications. 
Subcarrier  optical  communication  using  one  form  of  periodic  orbits  commonly  known  as 
self- sustained  pulsations  (SSP)  [5]  has  been  experimentally  demonstrated  both  for 
analog  [6]  and  digital  [7]  applications.  The  signal-to-noise  ratio  (bit-error-rate)  of  the 
analog  (digital)  subcarrier  systems  depends  on  the  nature  and  magnitude  of  the  noise  of 
the  subcarrier.  Since  only  one  harmonic  in  SSP  is  used  as  the  subcarrier,  the  amplitude 
and  phase  noise  of  each  harmonic  in  SSP  needs  to  be  characterized.  In  this  repoty  we 
report  the  first,  to  the  best  of  our  knowledge,  calibrated  measurement  of  amplitude  and 
phase  noise  of  self-pulsations  in  two-section  Fabry-Perot  laser  diodes  (LDs). 

4.2  Experimental  Setup 

Figure  1  illustrates  the  experimental  setup  for  the  measurement  of  amplitude  and  phase 
noise  of  SSP  in  commercially  available  780  nm  SSP  LDs.  The  SSP  is  first  converted  into 
RF  signals  with  a  high-speed  photodetector  (3  dB  bandwidth  >14  GHz).  The  RF 
spectrum  of  SSP  is  monitored  using  a  broadband  RF  spectrum  analyzer  (HP  8565E).  By 
using  bandpass  filters,  individual  harmonic  of  the  RF  signal  can  be  selected  for 
amplitude  and  phase  noise  measurement  using  the  HP  3048A  Noise  Measurement 
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System,  which  is  a  calibrated  measurement  system  dedicated  for  RF/microwave  oscillator 
noise  characterization  for  frequency  offsets  as  small  as  1  Hz. 


Figure.  5  Experimental  setup  for  the  amplitude  and  phase  noise  measurement 


4.3  Results 


RF  Spectrum 
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Figure.  6  RF  spectrum  of  SSP  with  a  1  GHz  fundamental  frequency  in  a  laser  diode. 
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A  typical  RF  spectrum  of  SSP  is  shown  in  figure  2.  At  70  mA  DC  bias  current,  the  SSP 
frequency  is  around  1  GHz  (The  threshold  is  at  50  mA).  More  than  15  harmonics  can 
be  observed.  Each  harmonic  has  a  certain  linewidth  which  increases  with  its  order. 
Details  around  the  fundamental  harmonic,  given  in  the  inset,  reveals  a  3-dB  linewidth  on 
the  order  of  1  MHz.  Because  the  fundamental  harmonic  was  used  in  previous 
experiments  [6,7]  and  because  of  the  instrumentation  frequency  limitations,  we  focused 
on  the  noise  characterization  of  the  fundamental  harmonic  in  the  present  effort.  With  a 
1  GHz  bandpass  filter  of  400  MHz  bandwidth,  the  amplitude  and  phase  noise  of 
fundamental  harmonic  are  investigated  nearly  independent  of  each. 


Figure.  7  Amplitude  noise  spectrum  of  the  fundamental  harmonic. 


The  amplitude  noise  is  measured  first  because  the  AM  noise  interface  offers  a 
better  than  110-dB  discrimination  against  phase  noise.  The  amplitude  noise  spectral 
density  in  dBc,  which  represents  the  AM  noise  sideband  to  carrier  power  ratio,  is  given  in 
figure  3.  Comparing  the  RF  and  AM  noise  spectra,  it  is  apparent  that  amplitude  noise 
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has  a  negligible  contribution  to  the  overall  noise  power  leading  to  the  1  MHz  linewidth. 
Therefore,  the  noise  of  the  fundamental  harmonic  has  to  be  primarily  phase  noise.  The 
phase  noise  interface  is  essentially  a  delay-line  frequency  discriminator  which  measures 
the  noise  frequency  deviation  from  the  center  frequency  due  to  phase  noise.  Figure  4 
shows  the  phase  noise  spectrum  of  the  fundamental  harmonic  for  1  Hz  to  40  MHz 
frequency  offset  from  the  center  frequency  of  the  carrier.  Below  1  kHz  offset,  the  noise 
frequency  deviation  is  1/f  in  nature.  The  peak  around  60  Hz  results  from  the  power 
supply.  For  the  frequency  offsets  above  10  kHz,  frequency  deviation  noise  is  white  with 
a  magnitude  of  around  600  Hz  /  yffiz  .  The  decrease  near  40  MHz  offsets  is  due  to  the 
bandwidth  limitation  of  the  HP  3048A  system.  Assuming  a  600  Hz  1 4Th  pure 
frequency  fluctuation  which  is  white  in  nature,  the  theoretically  [8]  estimated  carrier 
linewidth  is  about  1.13  MHz  which  matches  very  well  with  the  observed  RF  spectrum. 


Frequency  offset  (Hz) 


Figure.  4  Phase  noise  spectrum  of  the  fundamental  harmonic 


4.4  Summary 


We  have  investigated  experimentally  the  noise  behavior  of  self-sustained  pulsation  in 
laser  diodes.  Frequency  deviation  and  amplitude  fluctuation  of  the  pulsation  have  been 
characterized  in  the  range  of  1  Hz  to  40  MHz.  Even  though  phase  noise  dominates  over 
amplitude  noise,  both  are  1/f  in  nature  for  offsets  below  1  kHz,  and  white  for  offsets 
above  10  kHz.  This  technique  can  be  applied  to  the  noise  measurement  of  the  high- 
frequency  oscillation  in  ECLDs  also..  In  addition  to  understanding  and  ultimately 
improving  the  performance  of  periodic-orbit  subcarrier  optical  communication  systems, 
further  research  along  this  line  will  be  useful  for  other  applications  such  as  all-optical 
clock  recovery  [9] . 

5.  Conclusions 

In  conclusion,  we  have  determine  the  underlying  mechanism  of  the  high-frequency 
intensity  oscillations  in  laser  diodes  with  short  external  cavity.  From  a  nonlinear  dynamic 
point  of  view,  Hopf  bifurcation  is  responsible  for  the  occurrence  of  the  high-frequency 
intensity  oscillations.  However,  Hopf  bifurcation  is  only  the  necessary  condition.  The 
necessary  and  sufficient  condition  are  that  Hopf  bifurcation  occurs  and  the  limiting  point 
of  existence  of  a  higher-order  mode  is  above  the  stability  boundary  of  a  lower-order 
mode.  The  noise  characteristics  of  the  high-frequency  oscillations  can  be  measured 
using  the  HP  3048A  phase  and  amplitude  noise  measurement  system. 
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Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 

a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Material 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence, 
reliability  science,  electro-magnetic  technology,  photonics,  signal 
processing,  and  computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


